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Abstract. We present a study of magnetic field induced quantum phase transitions 
in insulating systems. A generalized scaling theory is used to obtain the temperature 
dependence of several physical quantities along the quantum critical trajectory {H = 
He, T — !■ 0) where is a longitudinal external magnetic field and He the critical 
value at which the transition occurs. We consider transitions from a spin liquid at a 
critical field Hci and from a fully polarized paramagnet, at Hc2, hito phases with long 
range order in the transverse components. The transitions at Hci and Hc2 can be 
viewed as Bose-Einstein condensations of magnons which however belong to different 
universality classes since they have different values of the dynamic critical exponent. 
Finally, we use that the magnetic susceptibility is an entanglement witness to discuss 
how this type of correlation sets in as the system approaches the quantum critical point 
along the critical trajectory, H = Hc2, T — > 0. 
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1. Introduction 

Recently there has been an intense study of field induced quantum phase transitions in 
metalhc [1112113111131^1 and insulating materials [711111311011111 ?]• These field induced 
transitions are generally associated with a Bose-Einstein condensation of magnons. In 
the insulating case, most of the experimental studies, have concentrated in obtaining the 
shift exponent that characterizes the shape of the critical line in the neighborhood of the 
quantum critical point (QCP). However, it is important to have additional experimental 
information that can be used to fully characterize the universality classes of the different 
zero temperature field induced transitions. For this purpose it is convenient to obtain the 
thermodynamic behavior of the system along a special trajectory in the phase diagram. 
This quantum critical trajectory consists essentially in sitting at the quantum critical 
point, by fixing the magnetic field at its critical value and vary the temperature. Field 
induced transitions are generally associated with soft modes. Then, at the QCP where 
the gap for excitation vanishes, physical quantities as specific heat, magnetization and 
susceptibility have power law temperature dependencies determined by the quantum 
critical exponents. Our aim here is to obtain this behavior. This strategy has been 
intensively explored in the study of heavy fermion materials, in this case, fixing the 
pressure at its critical value for the disappearance of magnetic order ^3]. This approach 
is particularly useful for insulating materials where the magnons that become soft at 
the QCP account for most of the low temperature thermodynamic behavior. In metals, 
even in field induced transitions, one has to consider also the contribution of charge 
carrying excitations [H 12113 11113101 • The soft magnon modes couple to these excitations 
and both may be strongly renormalized. 

We will consider two types of field induced transitions. Firstly, we study the 
zero temperature transition in an antiferromagnet system, from the fully polarized 
paramagnetic state into a phase with transverse components of the magnetization as the 
longitudinal magnetic field is reduced to a critical value Hc2 [El ■ This is a second order 
transition which is characterized by a dynamic critical exponent z = 2 It can be 
identified with a Bose-Einstein condensation of magnons with the magnetic field playing 
the role of the chemical potential. We also consider the transition from a disordered 
spin-liquid phase to an antiferromagnetic phase ^T], at a critical magnetic field, Hci- 
In this case the field is increased to the critical value Hqi at which the the singlet 
gap vanishes and magnetic long range order sets in. This transition is characterized 
by a dynamic critical exponent z = 1 ^3]. In both cases the transition is approached 
from the disordered side, i.e., with H > Hc2 and H < Hqi- Although these transitions 
have been intensively studied [3 HDl HH 1121 HE] the physical behavior along the quantum 
critical trajectory, H = He, T — »• has not been sufficiently characterized. In this paper 
we fully determine the thermodynamic properties and the appearance of entanglement 
along this line. We hope in this way to motivate further experiments that go beyond 
obtaining the shift exponent of the critical line. Notice that away from the QCP, in 
the disordered phase, the excitations are gapped and the physical behavior is thermally 
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activated. 



Our results are obtained using a generalized quantum scaling theory which was 
previously applied to heavy fermion materials ^3 ^] . This theory describes completely 
the thermodynamic behavior along the quantum critical trajectory. The reason is that 
we consider three dimensional materials and the effective dimensions fri\ IT^ at the field 
induced quantum phase transitions dejf = d + z are at, or above, the upper critical 
dimension = 4 for these transitions. In this case, the quantum critical exponents 
assume Gaussian or mean-field values, being fully determined. Logarithmic corrections 
which arise in the case dg// = dc are not obtained in the present approach. The 
microscopic justification of the scaling theory is provided by renormalization group 
calculations P^ ll6j. In the last section we discuss how entanglement sets in among the 
spins along the quantum critical trajectory. 

2. Scaling analysis of the quantum phase transition 

We consider first an antiferromagnetic system with longitudinal anisotropy and a strong 
magnetic field applied along the easy axis direction. We study the T = transition, 
with decreasing magnetic field, from the fully polarized paramagnet to a phase with 
long range order in the transverse components of the magnetization at a critical field 
Hc2{T = 0) P3]. This problem can be treated as a dilute gas of bosons with the effective 
action given by , 



where 6 = H — Hc2- The field ipO^y^) is a two-component field representing the 
components of the spins transverse to the direction of the magnetic field and vq takes 
into account the spin-wave interactions. The transition at 5 = has a dynamic exponent 
z = 2 due to the ferromagnetic-like dispersion of the magnons [14J. This action has been 
extensively studied in the context of the non-ideal Bose gas [201 , the dilute Bose gas 
PT] and superfiuid- insulator transitions In our case it is useful due to the small 
number of magnons excited at the field induced quantum phase transition j2H] . It takes 
into account the ferromagnetic-like dispersion of these modes ^3] and incorporates a 
constraint in their total number [23j . 

The thermodynamic properties of systems, close to the QCP described by the action 
in Equation (^, can be obtained from the free energy density. This has the scaling form 



near the QCP. The zero temperature critical exponents a, u and the dynamic exponent 
z are related to the dimensionality of the system d by the quantum hyperscaling 
relation, 2 — a = v{d + z) jTHj. In general for dejf = d + z > 4, i.e., above the 
upper critical dimension dc = 4, the exponents associated with the QCP at 5 = 




(1) 




(2) 
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take Gaussian values and in particular the correlation length exponent, u = 1/2. For 
the action in Equation these exponents remains Gaussian even below dc = 2 |2Uj . 
Although within the renormalization group approach for d = 3, the transition at 
(5 = is controlled by the Gaussian fixed point, the spin- wave coupling vq acts as a 
dangerously irrelevant interaction for deff > 4 and must be dealt with carefully |19j . 
Perturbation theory in powers of Vq leads to a temperature dependent critical line given 
by, 6{T) = Hc2{T) - Hc2 + VoT^''^ = with the shift exponent ip = z/{d + z - 2) = 2/3 
in three dimensions pTIj . 

In order to obtain the correct scaling behavior near the QCP, the scaling function 
F{x) in Eq. El must have the asymptotic behaviors, 

{constant for a; ^ 
(3) 
for X — s> oo 

The first guarantees that we recover the correct behavior at T = 0, with / oc 
The second with p = {a — a)/ uz, yields 

fiT) ^ AiTMT)\'-^ (4) 

where A{T) = j24j. Tilde exponents refer to finite temperature transitions and a 
is the Gaussian thermal specific heat exponent of the 3d — XY model. This is related 
to the thermal Gaussian correlation length exponent v through the Josephson relation, 
2 — a = vd. Since v = i) = 1/2 for c? + z > 4, the free energy in the neighborhood 
of the QCP can be written as, / oc T\5{T)\'^'^ {p = 1). From this expression we obtain 
the temperature dependence of the magnetization, susceptibility and specific heat at 
H = Hc2, 

m =aiT'/^ + hv'J^T'/^ (5a) 
X =a,T'/' + h,{l/v,f/'T'l' (55) 
Cy = a,T^'^ + hvl"T'" (5c) 

respectively, where the {ai,bi) are constants. For each quantity, the first term is the 
Gaussian contribution and the second arises from corrections due to the dangerous 
irrelevant spin- wave interaction vq. Notice that, except for the susceptibility and 
correlation length (see Table Q), for which vq behaves as a truly dangerously irrelevant 
interaction as it appears in a denominator, the purely Gaussian term is dominant at 
low temperatures. A similar analysis can be carried out for the transition from the 
spin liquid to the antiferromagnet at Hci with a dynamic exponent z = 1. The results 
are given in Tabled This case is marginal, since d^ff = dc = 4: is the upper critical 
dimension and there are logarithmic corrections to the quantities in Table Q These 
however are not obtained in the scaling approach (see Refs. I25j). 

Away from the QCP, for H > Hc2 and H < Hci, there are gaps for excitation 
of spin- waves and the crossover temperature Tx oc l^l*^^ gives the energy scale for 
the thermally activated behavior of the thermodynamic functions in these regions of 
the phase diagram. Since these gaps vanish as 161"^ their measurement allows for a 
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Physical Quantity 




Hci{z = 1) 


Hc2{z = 2) 


Shift exponent 




1/2 


2/3 


Magnetization 


-Of /OH 




2^3/2 


Susceptibihty 


-d^f/dH^ 


1 


1 yl/4 


Specific heat 


-Td^f/dT^ 




T3/2 


Correlation length 






1 y-3/4 



Table 1. Power law temperature dependence of physical quantities along the critical 
trajectory [H = He, T ^ 0). Logarithmic corrections are not included. Only the 
dominant contribution in the low T limit is given |26j. 



determination of the gap exponent uz. Notice that the difference in universality classes of 
the two transitions studied above arises essentially from the different dispersion relations 
of the soft magnons at the QCP. 

Finally at T = in the ordered phase, it is also necessary to take into account the 
dangerous irrelevant spin- wave interaction vq. The scaling form of the free energy is 
[IH], / oc \S\'^^'^~^^^ [uol^l*-'^^^"^-*^^] . The dangerous irrelevant nature of vq is manifested 
in the fact that the scaling function Fy[x —>■ 0] o(. 1/x. This yields 

J ^ ^^|5|(d+.-4)/2 - — 

for d+z > 4. This mean-field behavior gives rise to a /on(/«tMdma/ magnetization varying 
linearly with the distance to the QCP, i.e., M oc \H — Hc\- 



3. Entanglement and the quantum phase transition 

Recently, there has been a large interest in characterizing entanglement in systems near 
quantum critical points and in macroscopic magnetic systems [27]. For a system with 
N spins, Wiesniak et al. [2H] have shown that the magnetic susceptibihty acts as an 
entanglement witness and that whenever, 

X = Xx + Xy + Xz< -jij; (6) 

there is entanglement between the individual spins of magnitude /. The Xi cire the 
susceptibilities along three orthogonal axis measured in the same quantum state. A 
quantum complementarity relation involving the susceptibility and magnetization can 
also be obtained ^| (T 7^ 0), 

kTx 

1 H < 1 (7) 

Nl - 

This is particularly useful when applied to low dimensional materials, or systems at 
quantum criticality, as temperature can be reduced without any phase transition. In 
the equation above, it is useful to define the quantity E{T,H) = 1 — (kTx)/Nl 
which provides a measurement of entanglement, while the last term S = M'^/N'^l'^ 
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represents local properties. Equation [7| shows the interplay between these quantities 
since, 0<E + S<1. 

It is interesting to apply the relations above to the previous problems. Let us 
consider the case z = 2, with the applied magnetic field fixed at the critical value Hc2- 
For T = and H > Hc2 the system is in a fully polarized state, M = Nl, and we must 
have, E{T —^0,H = Hc2) — > 0. This implies that, as T — > 0, x = Xx+Xy+Xz = Nl/kT. 
Since the system is already fully polarized at T = 0, we must have XziT H = 

Hc2) and consequently, Xx{T ^0,H = Hc2) = Xy{T ^0,H = Hc2) = Nl/2kT. 
Assuming that the transverse uniform susceptibilities have already taken their low 
temperature asymptotic behavior as entanglement sets in with decreasing temperature, 
an approximate condition for the appearance of this type of correlation can be obtained 
from Equation ^ in terms of the longitudinal susceptibility alone. This is given by, 
Xz < Nl/2kT. For the specific case that the spin / = 1, this condition can be made 
precise [211. It turns out that whenever, Xz(T, H = Hc2) < {9/16){N/kT) entanglement 
sets in among the spins. 

As temperature increases along the line H = Hc2 the entanglement measure E{T) 
also increases. Since the magnetization for H = Hc2 decreases as M = A^/(l — aT^/^), 
the complementarity relation Equation ((Tj) implies that, 

E{T) = 1-^<1~{1- aT'/^f. (8) 

This can be written as, E{T) < f{T/Tc) where Tq = (l/a)^"^^. The quantity a is easily 
calculated and we get (/ = 1), 

^ ^^^^ ^ 

^ ~ C(3/2)2/3 mkB t;2/3- ^ ' 

This is just the Bose-Einstein condensation temperature of a system of bosons. In 
this expression, m = {h^ /2D) is the mass of the magnons with spin-wave stiffness D 
and V = V/N, with A^ the total number of spins in the volume V. In our problem 
described by Equation this arises due to the constraint in the number of modes 
and that the dominant contribution for the decay of the magnetization is the Gaussian 
one, the interaction vo being truly irrelevant for this quantity (see Table Q). Thus at 
low temperatures where the spin-wave approximation holds entanglement scales with 
the characteristic temperature Tc at the quantum critical point. This is an interesting 
feature of this quantum phase transition associated with a soft mode. Although the 
crossover temperature Tx oc \6\'^^ vanishes at the QCP and excitations become gapless, 
the mass of the bosons remain finite providing an energy scale even at the QCP. 

It is useful to explore the analogy of the present magnetic problem with a true Bose- 
Einstein condensation of bosonic particles [30! to gain insight in the latter problem. 
The relevant Bose-Einstein transition in this case is the zero temperature density- 
driven transition in a system of interacting bosons from the incompressible insulating 
phase to the superfiuid [221 I2H1 IHI]- The control parameter 6 is given hy ji — nc 
where iic is the critical (interaction dependent) value of the chemical potential. The 
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magnetization in the magnetic problem corresponds to the number of condensed bosons 
and the longitudinal susceptibility to the compressibility defined as k = —d'^f/dfi^. 
The transverse uniform susceptibility can now be associated with the order parameter 
susceptibility of the superfluid [301 and diverges for T — at the QCP |22], A* = f^c- The 
analogy with the magnetic case yields a criterion for the appearance of entanglement 
along the quantum critical trajectory (/i = fic, T — * 0) that can be expressed solely 
in terms of the compressibility. For (/ = 1) this is given by, k < {9/16){N/kT). The 
characteristic temperature which constraints the entanglement measure is the Bose- 
Einstein critical temperature of bosons with density n{^c)- Entanglement in this case 
implies the establishment of phase coherence among the particles jSHj. 

4. Conclusions 

We have obtained the thermodynamic properties at the quantum critical point of 
magnetic field induced phase transitions using a scahng approach. We considered first 
the zero temperature transition from the saturated paramagnetic phase to a phase with 
long range order in the transverse components of the magnetization with decreasing 
field. We also presented results for the field induced transition from the spin liquid 
to the antiferromagnet. These transitions are in different universality classes, due to 
the distinct values of the dynamic critical exponents. These are determined by the 
dispersion relation of the gapless excitations at the QCP. Since for 3d systems and the 
dynamic exponents considered here deff > dc, the critical exponents associated with 
both QCP can be immediately obtained in this case. Although for deff > dc the fixed 
point governing the quantum phase transition is Gaussian, in both cases, the quartic 
term in the action due to spin-wave interaction is dangerously irrelevant and must be 
considered. It plays a fundamental role in determining the critical line and the behavior 
along the quantum critical trajectory of the correlation length and susceptibility. 

For the problem described by the action Equation (^, the zero temperature 
exponents remain Gaussian even below d^ = 2 [20j. However, as discussed below 
Equation (j3)) the behavior along the quantum critical trajectory in the presence 
of dangerously irrelevant interactions it is also affected by the thermal exponents. 
For this reason we restricted our analysis to 3d systems. Finally, we have 
investigated entanglement properties along the critical line as witnessed by the magnetic 
susceptibility. We have shown that at the QCP, the Bose-Einstein condensation 
temperature provides a well defined energy scale for the appearance of entanglement 
among the spins 
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